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^' Abstract 

^ . 
JL , We combine earlier investigations of linear systems with Levy fluc- 

C^ I tuations [Physica 113A, 203, (1982)] with recent discussions of Levy 

J^ ■ flights in external force fields [Phys.Rev. E 59,2736, (1999)]. We give 

a complete construction of the Ornstein-Uhlenbeck-Cauchy process as 
a fully computable model of an anomalous transport and a paradigm 
^ ' example of Doob's stable noise-supported Ornstein-Uhlenbeck pro- 

» . cess. Despite the nonexistence of all moments, we determine local 

characteristics (forward drift) of the process, generators of forward 
and backward dynamics, relevant (pseudodifferential) evolution equa- 
tions. Finally we prove that this random dynamics is not only mixing 
(hence ergodic) but also exact. The induced nonstationary spatial 
process is proved to be Markovian and quite apart from its inherent 
discontinuity defines an associated velocity process in a probabilistic 
sense. 



1 Motivation 

The casual understanding of the central limit theorem (in reference to the 
Boltzmann-Gibbs thermostatics), combined with the need to have clearly 
specified the mean features (moments and local conservation laws) of ran- 
domly implemented transport, at or off thermal equlibrium, resulted in an 
obvious predominance of Gaussian laws of probability and diffusion processes 
in typical statistical (eventually probabilistic, cf. the omnipresence of the 
Brownian motion conceptual background) analysis of physical phenomena. 

Presently, we observe a continually growing recognition of the profound 
role (ubiquity, [![]) played by non-Gaussian Levy distributions (probability 
laws) in both a consistent probabilistic interpretation of various experimen- 
tal data and in a stochastic modelling of physical phenomena, followed by 
numerical and realistic experimentation attempts to verify (or rather falsify) 
probabilistic hypotheses. 

Generically, Levy's probability laws appear in the context of anomalous 
diffusions (mostly subdiffusions that are modelled in terms of continuous ran- 
dom walks, 0). On the other hand, under the name of Levy flights, 0, Q we 
encounter stochastic jump-type processes which are explicitly associated with 
those distributions. That allows in turn to model quite a variety of trans- 
port processes, cf. which are either regarded as (non)typical phenomena 
of nonequilibrium statistical physics or as manifestations of a complex non- 
linear dynamics with signatures of chaos, yielding an enhanced diffusion in 
particular. 

We focus our attention on Levy flights which are considered as possible 
models of primordial noise, 0, Q. (Wiener noise or process is normally in- 
terpreted to represent a statistical "state of rest" of the random medium). 
Generically, the variance and higher cumulants of those processes are infi- 
nite (nonexistent). There is also physically more singular subclass of such 
processes for which even the first moment (mean value) is nonexistent. Thus 
we need to relax the limitations of the standard Gaussian paradigm: we face 
here a fundamental problem of establishing other means (than variances and 
mean values) to characterise statistical properties of Levy processes (frac- 



tional moments of Ref. |^ are insufficient tools in this respect). 

Specifically, if a habitual statistical analysis is performed on any exper- 
imentally available set of frequency data, there is no obvious method to 
extract a reliable information about tendencies (local mean values) of the 
random dynamics. Nonexistence of mean values and higher moments may 
also be interpreted as the nonexistence of observable (e.g. mean, like drifts 
or local currents) regularities of the dynamics. Moreover, the jump-type pro- 
cesses usually admit arbitrarily small jumps (with no lower bound) and finite, 
but arbitrarily large jump sizes (with no upper bound). Any computer simu- 
lation utilizes both the lower (coarse-graining) and upper bound on the jump 
size, p, ^, |10|, and any experimental data collection involves such limitations 
as well. Mathematically, that puts us in the framework of standard jump 
processes for which the central limit theorem is known to hold true in its 
Gaussian version (even if we account for the abnormally slow convergence to 
a Gaussian, in view of long tails of the probability distribution, |^). There- 
fore, there is no clear-cut procedures allowing to attribute an unambigous 
statistical interpretation in terms of Levy processes to given phenomenolog- 
ical data. In a drastic contrast to a traditional Gaussian modelling. Mere 
scaling arguments, reflecting the self-similiar patterns of sample paths, are 
insufficient as well. 

Although no realistic formulation of a fluctuation-dissipation theorem is 
possible in that case (nonexistence of variances), we can give a meaning 
to a theory of Levy flights in external force fields, 0, under a simplifying 
assumption that force fields define linear processes with Levy fluctuations. 
The corresponding velocity processes were introduced in Ref. |0 ( see also 
irrp, but we shall give a complete construction of the related jump- type 



stochastic process, together with a detailed characterization of the dynamics 
of induced spatial displacements. Our strategy is thus substantially different 
from that typically followed in the current literature, 0. For example, the 
configuration space Langevin equation, 

^ ^ ^ + rm (1) 

at m7 

where m is the mass of transported particle, 7 stands for the friction con- 



stant and r] represents any conceivable generalization of the white-noise that 
employs Levy stable statistics, 0, corroborates a tacit assumption that some 
kind of the the standard Smoluchowski projection (the large friction limit, 
normally employed in the Brownian motion context, |]12|) from the phase- 
space to spatial only dynamics is possible. This is certainly not realizable in 
the non-Gaussian Levy case. 

Another delicate question is to settle possible physical origins of the spa- 
tial noise. That issue seems to be conceptually easier to handle on the veloc- 
ity/momentum space level. However, another delicate problem is immediate: 
in case of the Brownian motion (Ornstein-Uhlenbeck process) spatial trajec- 
tories were by construction differentiable to give meaning to the velocity 
concept (even though accelerations were nonexistent anyway, somewhat con- 
flicting with the naive but widely spread usage of the white noise-supported 
Langevin equation as the second Newton law analogue). The conseuqent 
exploitation of Levy processes with their intrinsic discontinuities, seems to 
set an unresolvable obstacle in this (velocity notion) respect. 

We shall demonstrate that this is not literally the case. For example, we 
can prove that the spatial random variable x(t) of the Ornstein-Uhlenbeck- 
Cauchy process cannot possess derivatives in the sense of standard mathe- 
matical analysis, nonetheless this process has derivatives in a weaker, proba- 
bilistic sense. Hence, it is legitimate to interpret u(t) as a velocity analogue 
attributed to the instantaneous (random) location x(t), though this notion is 
more distant from classical intuitions than the velocity variable of the stan- 
dard Ornstein-Uhlenbeck process (not differentiable, hence not yielding any 
analogue of a Newtonian acceleration). 

Our analysis departs from a generalization of the Ornstein-Uhlenbeck 
process due to Doob, [0, where a symmetric Levy (stable) noise was assumed 
to take place of the standard Wiener noise. A complete description of a 
concrete, computable in full detail Ornstein-Uhlenbeck-Cauchy process (with 
a familiar lorentzian as a probability law for velocity displacements) is our 
principal goal in the present paper. In addition we shall pay attention to 
intrinsic complications of the random dynamics by investigating a standard 
chain of its possible features (ordered with respect to the complication level): 



ergodicity, mixing and exactness. 



Langevin equation with a linar (harmonic) 
force and Cauchy noise 



The starting point for Ornstein and Uhlenbeck,||13|, |T4|, was the dissipative 
Langevin equation 

^ = -Au(t) + A{t) (2) 

where u(t) is a random variable describing the velocity of a particle, A > 
is a dissipation constant, and A{t) is another random variable whose proba- 
bilistic features are determined by the probability distribution of u(t), which 
is assumed to satisfy a concrete law when t ^ oo. Because u(t) may have 
no time derivative, equation (2) was soon replaced by another one, namely 

du{t) = -Xu{t)dt + dB{t), u(0) = uq (3) 

which received a rigorous interpretation within the framework of stochastic 
analysis |0. In the case when the probability distribution of u(t), t -^ oo, 
is the Maxwell one, B{t) must be a Gaussian process, and the formula (3) 
leads to the classical Ornstein-Uhlenbeck process. 

Here, we discuss properties of the process u(t), and the corresponding 
process of displacements x(t), in the case when B = {B(t))t>o is the Cauchy 
process, that is when B satisfies the following conditions: 

a) B has independent increments, i.e. given ti < ... < tn, the differences 
B(t2) — B(ti), Bits) — B{t2),..., B{tn) — B{tn-i) are mutually independent 
random variables, 

b) B has stationary increments, i.e. the probability distribution of B{t + 
r) — B{t) is independent of r, 

c) B is continuous in probability, that is limits B{t) = B{s) in probability, 

d) the characteristic function of B is given by 

where ipip) = <^^b|- 



All the above requirements form a mathematically consistent definition of 
the Markovian jump-type process in question, e.g. Cauchy process. Notice 
that a suitable modification of the condition d) (set ip{p) = —a^p^ in the 
exponent; we refer to the general form of the Levy-Khintchine formula) would 
leave us with the familiar Wiener process. 

From a physical point of view, solutions of induced partial (here, pseudo-) 
differential euqations are most important, and those incorporate transition 
probability densities and densities of the process. 

Notice that the process of displacements is determined by u(t) in the 
standard way 

t 
x(t) = x(0) + f u{T)dT, x(0) = xo (4) 



Hence, we should be able to derive relevant densities and transition densities 
not only for the velocity process but also for the induced spatial process. 
Additionally, if we wish to interpret u(t) as a genuine velocity field for the 
process of spatial displacements x(t) (the mere formal attribution of the 
velocity name to our random variable u(t) is highly misleading, in view of an 
apparent discontinuity of sample paths), a careful analysis of differentiability 
properties (in what sense ?) of x(t) is here necessary. 
By integrating equation (2) we obtain that for t > s 

t 



u(t) = e-^(*-")u(s) + e"^* / e^"rf5(r) (5) 

s 

The integration of (3) yields 

t T 

x(t) = x(s) + f[e'^^^^'^u{s) + e-^^ f e^'^dB{f3)]dT 

s s 

= x(s) + u(s) - J dr{^' J e^^dBiP) 

s s 

Integrating the last summand by parts and using the double integration 
formula involving a derivative with respect to the interior integral, cf. (3.12) 



in 0, we get 

t 



t ^ T \j. t t 



1 _ e-A(i-s) }l - e"^(*"^) 



Hence 

x(t) = x(s) + ^ ^ u(s) + J dB{T) (6) 

s 

which mimics (is identical with respect to the form) a fairly traditional ex- 
pression for a spatial random variable of the standard Ornstein-Uhlenbeck 
process (with Wiener increments put instead of the Cauchy increments in 
the last summand). 

3 Probability densities and transition proba- 
bility densities for u{t) and x(t) 

There is a number of (equivalent) procedures to deduce a probability density 

of the process u(t) from the Cauchy increments statistics, see [^ ^. We shall 

follow a direct probabilistic route. 

In order to simplify the notation we write P[X = x] for the density of 

the probability distribution of a random variable X, that is P[X G F] = 

/r P[X = x]dx for T C R. 

Suppose that / is a continuously different iable function such that /(r) > 
t 
0, and let X = / f{T)dB{T). Terms of this functional form are encountered 

s 

in formulas (4) - (6). 

The random variable X is the limit of the following sum 

E/(rfc)[5(r,+i) -Bin)] (7) 

fc=0 

where s = tq < ri < ... < r„ = t is the partition of the interval [s, t]. 

Because the process B has independent increments, the probability density 

of (7) is the convolution of densities of its summands which, since the process 

has stationary increments, are equal to 

1 (T^(rfc+i - Tfc) ^ 1 a^f{Tk)ATk 

vr/(r,)(^)2 + or4(r,+i - n^ ~ vrx^ + (a2/(rfe)Arfe)2 



Because the Fourier transform maps the convolution to multiphcation 
and 

we get 

n-l 

P[ij:firk)[Bin^,) - B{n)]) = x] 

k=0 
n—l n— 1 

= iYle^'^'W^^^^^^-rix) = (exp(-a2|p|5]/(r,)Ar,))^(x) 

k=0 k=0 

where /^ and /^ denote the Fourier transform and its inverse respectively. 
By taking the limit n ^ oo we obtain that 

t 
P[X = x] = {eM-cr'\p\ J f{r)dr)nx) 

s 

and so the general formula 

1 <y'}f{r)dr 

P[X = x] = ^— (8) 

^a;2 + (a2//(r)rfr)2 

s 

is valid. We shall exploit Eq. (8) repeatedly in below. 

Remark 1: An apparent generalization of the previous observation is 
posssible. Assume that B{t) is a Levy stable process with the characteristic 
function 

^(p, t) = exp(-(T^t|p|"), < a < 2 

Then, there holds 

t t 

P[j i\r)dB{T) = x] = {exp[-a'{J nT)dr)\pnnx) 



Presently we shall use the formula (8) to calculate transtion probability 
densities of processes u(t) and x(t). 



Let /(r) = e ^^^ '^\ Then, by equation (5), 

Put = uu(s) = v] = ,,^ ^ ,,^ ^ ,, 9 

where cr^(t — s) = ^(1 — e"'^^*"''^), see e.g. also 0. 

Since u(0) = Uq, the probabihty density of u(t) is given by 

P[n(t) = u] = ^4^^ — - (10) 

We now turn to the process x(t). Since u(s) is independent of B{t) for 
all t > s 0, it is also independent of the integral /^ f{T)dB{T). Therefore, 
the probability distribution of the sum 

1 _ p-A(t-s) ^1 _ -A(t-T) 

uis) + / ^ dBir) 



A ^ ' 7 A 

is the convolution of its ingredients 

-A 

T 



Let /(r) = - — ^^j — —. Because of 



r 1 

j f{r)dr = ^(e-^(*-^) - 1 + X{t - s)) 



by formula (8) there holds 



t 



1 _ g-A(t-r) 1 (£)2(e-A(t-s) _ ^ ^ ^(^ _ ^)) 



^7 A "^^^^^ ""' 7rx2 + (£)4(e-A{t-.) _ 1 + A(t - s))2 

(11) 

On the other hand, by (10), we have 



A ' ' ' tt{u - uoe~^'a{t - s)y + a^{s)a'^{t - s) 

(12) 
where a(t — s) = ^--^ — -• 

The Fourier transform of (11) and (12) are equal to 

exp[-(^)2(e-^(*-^) _ 1 + A(t - s))\p\] (13) 

and 

exp[— a^(s)a(t — s)\p\]exp[—iuoe~^''^a{t — s)p] (14) 

9 



respectively. 
Because of 

A A 

the mulitiplication of transforms (13), and (14) followed by taking the inverse 
Fourier transform of the result, gives us a transition probability density of 
the spatial process 

P[x(t) = ,|x(s) = x]= piy, t\x, s) = ^ (^_,_ J^;/,)). ^^.(,^,) 

(15) 
where 

git, s) = ijne-'' - e'^^ + A(t - s)) 



p—^s p~^t 



and 

h{t, s) ^ 

Finally, because x(0) = Xq, the probability density of the process x(t) is 
given by 



1 im^-'' - 1 + At) 

vr (x - xo - uo^^f^y + (f )4(e-^* - 1 + \ty 



Pi^it) =x] = rA^^.o . ,..., ' —TTT^ (16) 



Compare e.g. the corresponding formula for the displacements of the stan- 
dard (Wiener noise-supported) Ornstein-Uhlenbeck process, [|12| . 



4 Properties of the process u{t) 

Considerations of the previous sections may leave us with an impression that 
a construction of the Ornstein-Uhlenbeck process supported by Cauchy noise 
is in fact complete. We have in hands not only Ito type stochastic differential 
equations that are amenable to direct computer simulations, [|10|, |^, but 
also explicit expresions for probability densities and transition probability 
densities for both processes: u(t) and x(t). In case of Markov processes such 
data are known to specify the process uniquely, |15|. 



However, some alarm bells need to switched on at this point. The stan- 
dard (stationary) Ornstein-Uhlenbeck velocity process is Markovian (in the 

10 



Gaussian case the Ornstein-Uhlenbeck process is the only continuous in prob- 
abihty stationary Markov process, ||16|), but the induced (integrated) spatial 
process is not Markovian. Using an exphcit expression for the transition 
probabihty density it is easy to verify that the Chapman-Kolmogorov iden- 
tity does not hold true. Therefore, Markov property is normally attributed 
to a two-component, phase-space version of the Ornstein-Uhlenbeck process, 
r^ . In case of the Ornstein-Uhlenbeck-Cauchy process the situation is some- 



what different. 

4.1 Markovianess and stochastic continuity 

First of all let us notice that u(t) is a time-homogeneous (but not station- 
ary) Markov process. Markov property is clear from the construction since 
Chapman-Kolmogorov identity can be verified by inspection and it is a clas- 
sic observation that nonnegative and normalized functions which obey the 
Chapman-Kolmogorov equation are necessarily Markovian transition proba- 
bilities. 

Since the probability density (10) of the process depends explicitly on 
time, our process u(t) is not stationary. 

Remark 2: That needs to be contrasted with the standard (Gaussian 
and stationary) Ornstein-Uhlenbeck process features where the transition 
probability density is time-homogeneous, while the density of the process 
does not depend on time at all. Indeed (we consider one spatial dimension 
and utilize dimensional units) the transition density 

pM^,s) ^ (V2.B{l-e.ri-2,(«-,)l})-V^.e.p(-2|p£^Mz2l^) 

with s < t, has an invariant density p(x) = ('j/2nD)^^^'^ ■ exp{—'jx^/2D). 
The drift of the process reads b{x) = —72; and p solves the Fokker-Planck 
(second Kolmogorov) equation dtP = DA^p — V xipp). 

Now, we shall demonstrate an important property (mentioned before in 
connection with the Ornstein-Uhlenbeck process) of the so-called stochastic 
continuity which is a necessary condition to give a stochastic process an 
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unambigous status, |T^, |T^ 0. Namely, we need to show that for any e > 
the following equation is satisfied 

limP[|u(t) - u(s)| > e] = (17) 

t — >s 

This equation is equivalent to 

lim / pt{u\v)du = (18) 

\u—v\>e 



Because of 



1 e + v{l - e"^*) e - v{l - e"^*) 



f / I X 7 . J- ^ e + -ull - e "M 

/ pt{u\v)au = 1 [arctan -— harctan 

J TV (J (t) 



\u—v\>e 

and remembering that o-^(t) = ^{^ — e""^*), the stochastic continuity prop- 
erty does follow. 

It is perhaps not useless to emphasize that in typical Gaussian process 
investigations, stochastic continuity of the process is a necessary (but still 
insufficient) condition for the process to have continuous sample paths. Hence 
it is always explictly mentioned in the context of diffusion pocesses, |jl5|] . The 
Cauchy noise-supported process is surely not diffusive and its trajectories are 
discontinuous (jump- type) paths, 0, 0). 

4.2 Local moments in the Cauchy case: forward drift 
— the existence issue 

The nonexistence of moments of the probability measure in case of the 
Cauchy process is another source of difficulties, since the standard local 
characteristics of the diffusion-type process like the drift and the diffusion 
function (or coefficient) seem to be excluded in the present case. 

However, for the considered Ornstein-Uhlenbeck-Cauchy process, the no- 
tion of the forward drift of the process proves to make sense (!). We shall 
first discuss the drift issue for the process u(t). 

Let us start with the following definition. Suppose p{y, t|x, s), t > s, is 
a Markov transition function and let Xt be the associated Markov process. 
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Guided by the analogy with diffusion processes we say that the process Xt 
has a drift (in fact, forward drift) if the following limit 

fe^i:^ J iy - ^)v{y,Axi s)dy (i9) 

\y~x\<e 

does not depend on the choice of e > 0. If so, then its value depending only 
on (x, s) we denote by 6(x, s) and call it the drift coefficient. 

Clearly, if p is homogeneous in time, then the drift coefficient depends 
only on the variable x. Let us emphasize that in the above definition we do 
not require the process Xt to have finite moments. 

We claim that the jump-type Markov process u(t) has a (forward) drift 
which reads h{v) = —Xv. 



At first we calculate the indefinite integral 

1 f. . cr^(t)du 

' [u — V)- 



Substituting z = u — f e~'^*, we rewrite that integral as 

oH^ r zdz t; ^^, _ r a^{t)dz 

71 J z^ + a^{t) n^ ' J z^ + a\t) 



a\t) 



V 



log{z^ + (T^{t)) + -(e^^* - l)arctan(- 



Hence 

I = 4^1og[(« - ve-'r + ^\t)] + -(e-^* - l)arctan[ ^~j;"'' ] 
and consequenly the limit 

lim-^-^(log[(t; + e - ve-^'Y + a\t)] - log[(t; - e - ve-^'f + a\t)]) 

t^O t ZTT 

1- ^ V .. V + e — fe~^* ,v — e — fe~^* , 
+ hm — (e — Ijfarctan -— — arctan -— ) 

13 



, f ,7r TT, 

= - A-(- + -) = -Xv 

TT / / 

exists and is e- independent. This is the forward drift of the process u(t) 
which proves a consistency of the derived transition probabihty density with 
the stochastic differential equation (5). 

To our knowledge, such consistency check has never been performed be- 
fore in discussions of Levy flights and anomalous diffusion processes. 

On the other hand, it is well known that for Markovian diffusion pocesses 
all local characteristics of motion (conditional expectation values that yield 
drifts and variances) are derivable from transition probability densities, sup- 
plemented (if needed) by the density of the process, cf. |jl8[. We have demon- 



strated that, in the non-Gaussian context, the nonexistence of moments does 
not necessarily imply the nonexistence of local characteristics (drifts) of the 
process. 

As a consequence, once a formal definition is adopted of a stochastic dif- 
ferential equation whose deterministic driving term (functionally unrestricted 
drift) is subject to perturbations by Levy flights, the process may still posess 
local characteristics (forward drift) that are in turn derivable by means of 
its transition density. Our derivation in the Cauchy noise case is limited to 
linear functions of random variables (linear systems, |P, ^. Possible gener- 
alizations to stochastic differential equations with driving terms represented 
by nonlinear and possibly time-dependent functions need to be carefully ex- 
amined. 

This is an uncomfortable situation, since a formal computer experimen- 
tation may not indicate any inconsistency of the formalism. Even worse, 
the uncommented visualization effectively may convey misleading or entirely 
wrong messages if uncritically accepted. (The rigorous existence theorems 
available in the mathematical literature pertain to linear systems as well, 



T^, pO[ and extend to perturbations by general Levy processes.) 
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4.3 Markov generators and Kolmogorov (Fokker- Planck 
type) equations 

Once densities and transition probability densities have been obtained from 
the first principles, we can invert the problem (that is a commonly shared 
viewpoint in the physics-oriented research) and ask for differential (evolu- 
tion) equations obeyed by them. The Fokker-Planck equation is an obvious 
example in case of Markovian diffusion processes, while various forms of the 
Master equation were adopted to extend the standard jump-processes (Pois- 
son or more gener ally-point processes) framework to more singular step or 
jump-type ones. 

In the case of unperturbed (free) Levy processes basically all interesting 
(covering stable laws of probability) evolution equations were classified by 
means of Fourier transform techniques, ||, ^, ^, ^ ^ . A disregarded point 
was that in case of Markov processes a single (Fokker-Planck or Master equa- 
tion - type) evolution equation does not characterise the process uniquely. 
Both forward and backward evolution equations need here to be involved, 
cf. fl^, |l^. Except for Refs. |^, |^ no attempt was made to investigate such 
equations for deterministically driven Levy systems. 

To elucidate that issue, we shall next consider the generator of a Markov 
transition function Pt{y\x) for the Cauchy-perturbed process (cf. the previous 
Section). 

Let us recall that it is defined by 

oo 

{Lf)ix) = lim^l J pM^)fiy)dy - fix)] (20) 

— oo 

where the domain of definition consists of all functions / G Co(R), whose 
limit on the right hand side in (19) exists uniformly with respect to the 
variable x. 

It is worth noting that when the transition function is stochastically con- 
tinuous (see the previous section), then the corresponding semigroup Tf in 
Co(R) defined by 

oo 

{TJ){x) = f Pt{y\x)f{y)dy (21) 
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is strongly continuous, and so its generator L is densely defined. 

In such a case we can also define an adjoint semigroup T^ acting on the 
space of (probability) densities -L^(R, dx), 



iT:p){u) = J ptiu\v)piv)dv (22) 

— oo 

Its generator we denote by L*. 

Arguments of the present section involve a little bit of a mathematical 
formalism to stay in conformity with the classic work of Feller and Dynkin 
on evolution equaitons for Markov processes, cf. |^, |T^ for references. 

Suppose L is the generator of the semigroup associated with the process 
u(t) and let L* be its adjoint. 

We wish to demonstrate that 

L = Lo + bV (23) 

and 

L* = Lo - V(6-) (24) 

where Lq is the generator of the Cauchy process B (we have used an explicit 
notation Lq = |V| in Refs. |Q, ^, see also 0) and b{v) = —\v. 

To this end, we first observe that for pt{u\v) given by the formula (9) the 
associated semigroup Tt maps Co(R) to Co(R). Since pt{u\v) is stochastically 
continuous, Tt is strongly continuous. 

Next, we calculate the Fourier transform of equation (20). 



oo oo 



(Lfrip) = liml[ / / ''J^l f{u)e-^^^dudv - f{p)] 



-oo — oo 



Substituting z = ve ^^ — u^ dv = e^^dz, we obtain that 

{Lfnp) = 



oo 



lim-[e^* / ^^- exp(-zze^V) ^2 1 ^4/^n / due^^{-iue^'p)f{u) - f{p)] 



-00 —00 



z^ + a\t) 
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= limi[e^*exp(-a2(t)e^*b|)/(e^V) - f{p)] 

= -ct'IpIHp) + A/(p) + Xpf'ip) 
By taking the inverse Fourier transform and using the identity 

oo oo 

iPfTiv) = ^ J e''^pf'{p)dp = -^ J {e^'"'p)'hp)dp 

— oo — oo 

= -f{v) - tv{pfy{v) = -f{v) - vf{v) 

we arrive at 

Lfiv) = Lofiv) - Xvf'iv) 

where Lq = — cr^|V|. 
Hence L = Lq + bV. 
Because of 

oo 

L*{p){u) = jim-[y ptiu\v)p{v)dv - p{u)] 

— oo 

so, by similar calculations as above, we obtain 

L*{p){u) = Lop{u) - V{b{u)p{u)) 
That ends the demonstration. 

As a consequence, since Tf : Co(R) — > C'o(R') is strongly continuous, 
almost all paths of the process u(t) are cadlag, that is they are continuous 
from the right and have finite left-hand limits (see chap. II sec. 4 in [0, vol.11). 
There follows also from Eqs. (23 ) and (24) that the transition probability 
function of the process u(t) satisfies the backward equation 

^^^ = Lop,{u\-){v) + b{v)VMn\v) (25) 

and the forward equation (the Fokker-Planck equation analogue) 

^Pl^ = LoPt{-\v){u) - V4b{u)pt{u\v)] (26) 
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for definitions see, for example, chap. 15 sec. 4 in [Q. 

It is trivial to check by inspection that the transition probabihty function 
of the (free) Cauchy process obeys both those equations with b{u) set equal 
identically 0. 

4.4 Asymptotics: ergodicity, mixing, exactness 

We have in hands an explicit expression for the density of the process u(t). 
One of efficient ways to investigate the complexity of the involved random 
dynamics is not necessarily via a direct recourse to sample paths, but rather 
via studying asymptotic properties of probability densities, cf . [^ . 

Let us consider the asjTiiptotic properties of the process u(t). By direct 
calculations we check that the density 



a^(oo) 



PoW 



ayX 



-nv? + (t'^(oo) -nv? + a^/X^ 



is stationary with respect to the dual semigroup T^ . Therefore, irrespectively 
of its initial probability distribution, P[u(oo) = u] = Pq{u). 

Remark 3: Normally, if we get a convergence of a density in the asymp- 
totic (t — > oo) regime to a unique density, we say about an asymptotic 
stability. In case when for every initial density we get a set spanned by a 
finite number of densities, we say about an asymptotic periodicity. We may 
also have a situation that every initial density is dispersed under the action 
of a Markov operator. That is related to the concept of sweeping. 

Our knowledge of the explicit formula of the transition probability den- 
sity for the semigroup T/ allows us to examine its ergodic properties in a 
more detailed way. For example, a point-wise convergence of the Cauchy 
process transition density to the stationary Cauchy probability density was 
ivestigated in p. 

Let us recall (see [^ for the definition and more details) that a Markov 



semigroup T^ is mixing, if for any density p, T*p tends to a stationary density 
in a weak sense, and exact, if this limit holds in the L^-norm. 



18 



Hence exactness is a stronger property and implies mixing, ergodicity of 
the dynamics being a straightforward consequence. 

(It mighth be worth noting that strong mixing properties of the standard 
Ornstein-Uhlenbeck velocity fields were discussed and visualized by computer 



simulations in |23|] . The point-wise convergence of the probability density of 
the process to its stationary limit was established in [p^ .) 

The dynamics induced by Cauchy noise (and other stable noises) shows 
higher level of complications and is not only mixing, but also exact. We shall 



provide a demonstration of this property in the spirit of Ref. [^ . 

In fact, what we claim is that T^ is exact (hence both mixing and ergodic). 
Since po is a stationary density we have to show that 

lim \\pt - polli = (27) 

t — >oo 

where pt(u) = J pt{u\v)p{v)dv and p{v) is an arbitrary initial density. To 
this end we need an auxiliary lemma which comprises the most technical 
segment of the paper. 

Lemma 

For t — > oo, ||pt(-|f) — polli ~^ uniformly in v on compact sets. 

Proof: We shall show that VA^ G N Ve > 3to > such that 

Vt > to Vt; G [-AT, N] \\pt{-\v) - Polli < e 
Let us begin from 



oo 



II rM II 1 f I ^ (^) ^^M I. 



-oo 



oo 



- n J \u - ve-^'f + a\t) 



-oo 
oo 



— OO 
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+ 



a\oc 




1 




— 


m2 


1 

+ a 


T7 T du = 

■4(00) 


o-2(oo) — o 


-2(t) 


TX 


J 'm2 

— oo 


+ a4 


it) 




+ 


a2(oo) 


oo 




\u' 


— 


(u 


- t;e"^*)2| 


An, 1 




TT 


J [{u 

-oo 


— 


ve^ 


\t\2 


+ 


(x^{t)][u^ + 


^^W] 




a' 


-(oo)[a^( 


'oo) - 


a' 


\t)] 


oo 
f 




dtt 







TT _7^ K + a\t)\[u^ + a4(oo)] 

The first summand, denoted by Ji, equals to e^^^ /{I — e"^^) and so is 
less than e/3 for all t > ti, provided ti = ^log(| + 1). 

For t > -2|-2 we have o"^(t) > |o"^(oo) and so the third summand, denoted 
by Is, can be estimated as follows 

/3 < ^^[.^(oo) + Ama\oo) - At)] / ^— ^^ 



^ ^ _io K+ (^)'P 



< ^^ KT^ = 8e 

TT (7° OO 



-At 



Therefore, for any t > t^ = max(^ log 2, ^ log ^) ther holds I3 < e/3. 
Finally, we estimate the second summand, denoted by /2. At first let us 
notice that 

o- (00) /■ \v\e \2u — ve \ 

du 



TX J \{u - i;e-^*)2 + aHt)\\u^ + a^ft)] 

—00 

_;^^o'2(oo) f \x\dx 



—00 
where x = 2u — ve~^^. 

Hence, for all t > ti = ^log(l + 2^:^) there holds \v\'^e-'^^^ < Aa^{t) 
for all V G [— A^, A^], and so 

x.o-2(oo) /■ \x\dx 1 e^^* 



TT 



-00 



a;4 + (4a4(t))2 a2(oo) (^ _ g-At)2 



Therefore, I2 < e/3 for all t > t2 = max(t4, t^), where t^ is deter- 
mined by 3e~'^*^ = o'^(oo)(l — e^^^^')'^e. Thus ||pi(-|f) — po||i < e for all 
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t > to = niax(ti, t2, ts) what ends proof of the lemma. □ 

Now, we are ready to adddress the exactness issue for T*. 
We observe that 



oo oo 



Pt - Polli = / I / Pt{u\v)p{v)dv - po{u)\du 



oo oo oo 



{pt{u\v) - po{u))p{v)dv\du < / \\pt{-\v) - po\\ip{v)dv 



-oo — oo — oo 

TV 



< / Ibt(-I^) - Po\\ip{v)dv + 2 / p{v)dv 

-N [-Af,Af] = 

where [— A^, A^]^ denotes the complement of the interval [— A^, A^]. 

Let us choose A^ such that the second integral in the above is less than 
e/4, and next, in conjunction with Lemma, we choose to such that for all 
t > to, there holds ||pt(-|f) - po||i < e/2 for all v e [-N, N]. 

Then 

AT 

J \\Pt{-\v) - po\\ipiv)dv < - 

-N 

and so \\pt — po\\i < e, which complets the exactness demonstration for T*. 

Remark 4: In our considerations of the exactness issue the stationary 
Ornstein-Uhlenbeck-Cauchy process has been employed. This process is a 
direct Levy stable analogue of the standard Gaussian Ornstein-Uhlenbeck 
process . There is however an important difference, |2^. In the Gaussian 



case all stationary Markov processes are Ornstein-Uhlenbeck (which in turn 
is unique). In particular, the standard process coincides with its reverse 
(fully anticipating) version. It turns out that in the Cauchy case there exist 
at least two different stationary Markov processes, since the reverse one does 
not coincide with the forward (nonanticipating) one. This means that the 
so-called statistical inversion of the Markovian dynamics (cf. a discussion 
and references in the closing section of Ref. ||26|), in case of stable Levy pro- 



cesses, makes a distiction between the "time arrow" direction. See e. g. also 
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time reversal and time adjointness problems encountered in Refs. 0, ^]. Pre- 
sumably all that derives from the exactness property. In modern attempts to 
devise a theoretical framework for nonequilibrium statistical physics, based 
on exploitation of clasically chaotic systems, most irregular and " most irre- 
versible" dynamical phenomena are known to be generated by exact systems. 

We could as well consider Eq. (1) with a nondissipative A < factor. 
However, then the asymptotic properties of the semigroup T^ change in an 
essential way. Indeed, since now the transition probability density can be 

written as 

1 m(el'l* - 1) 

Pt[u\v) = -—^ 5 — 

^ ' ^ vr(M - veW^y + (g(el^l* - l))^ 
then for any A^ G N and any density p there holds 

N 

lim / pt{u\v)p(v)dv = 

-N 

In consequence T^ is sweeping, [|^. It means that T^ has no stationary 
density and, in consequence, there is no probability law at all for the t ^ oo 
limit of the process u(t). 



5 Properties of the process x(t) 

5.1 Markovianess and forward drift 

In case of the classic Ornstein-Uhlenbeck process, it is well known that the 
spatial random variable does not represent a Markov process. 

A little bit surprisingly, in the present (Cauchy noise) case, we can prove 
that x(t) is a Markov process which is (like the previous u(t)) stochastically 
continuous. Moreover, while being a discontinuous process, nonetheless it 
has a forward drift equal b{s) = UQe~^^. 

To show Markov property it suffices to check the Chapman-Kolmogorov 
equation for the transition function given by equation (15). That immedi- 
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ately follows due to the additivity properties of functions g(t, s) and h{t, s) 
which enter the formula for p{y, t\x, s) 

g{t, t') + g{t', s) = git, s) 

h{t, f) + h(t\ s) = hit, s) 

The stochastic continuity can be shown by a direct verification of the 
formula (18) which, in the nonstationary case, reads 

lim / p{y, t\x, s)dy = 

t >S J 

\y—x\>e 

Also, by direct calculations we check that the limit 

lim _ [ J (y - x)p{y, t\x, s)dy] 

\y-x\>e 

does not depend on the chosen e cutoff, and equals ■uoc"'^*. 

Remark 5: It is worth pointing out that Markov property of the pure 
spatial process x(t) is a distinguishing feature of the Cauchy noise. It does 
not hold for other a-stable Levy processes, in particular for the Gaussian one 
(the standard Ornstein-Uhlenbeck process). The reason of this exception is 
rooted in a particularly simple form of the probability distribution of the 
process /, f{T)dB{T) when a = 1, see e.g. our Remark 1. 

In contrast to the velocity process our spatial process is no longer time- 
homogeneous. In the inhomogeneous case instead of a one-parameter semi- 
group we have a two-parameter family of operators Tf ., defined by 

oo 

iTtJ){x) = f p{y,t\x,s)f{y)dy (28) 



which satisfy the composition rule Tt^t'Tf^s = Tt^s- 

Therefore, we can introduce a time dependent generator by the following 
formula 

oo 

(M(s)/)(a:) = lim- [ / piy, t\x, s)f\y)dy - f\x)] (29) 

— oo 
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In analogy with our previous considerations, we can readily identify an 
explicit form of the generator M(s). Namely, let us assume that p{y, t\x, s) 
is the transition function of the process x(t) and Tt^s are operators associated 
with this function. Then 

M{s) = -a\s)\V\ + b{s)V (30) 

where o"^(s) is as (8) and b{s) = uoe~^^ . 

Because the major steps of the demonstration are essentially the same as 
in the case of the process u(t), we skip them here. 

Furthermore, let us notice that in view of lims_^oo ^(s) = xl^l ^*-'' 
for large t (i. e. asymptotically) the process x(t) converges to the Cauchy 
process with the transition function given by 

1 ta^/\ 

Pt{y\x) - 



n{y - a;)2 + t^a^X"^ 



Hence, the dissipation constant A > does the job "as usual", though 
with no recourse to the standard Maxwell-Boltzmann notion of thermal equi- 
librium and fluctuation-dissipation theorems. 



5.2 Sample paths features 

Let us turn to a brief discussion of the properties of sample paths of the 
process x(t). (We remember that sample paths of the jump- type process 
u(t) were cadlag.) 

In his seminal paper, 0, Doob proved that the displacements of the 
standard Ornstein-Uhlenbeck process satisfy 

hm "^^(^) - "^^(0) = (31) 

almost surely, i.e. the above limit holds for almost all sample paths (with 
probability 1). That is interpreted as an ergodic theorem applied to the 
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velocity process to give the strong law of large numbers, , and at the same 
time as a statement about the (sample) path of a single particle. 

In the final remark on p. 369, he also concluded that Eq. (31) holds true 
also in the case when the noise 5 is a stable process with the characteristic 
yU > 1. Hence, for the Cauchy process as well. 

However, this conjecture appears to be wrong, in view of the estimates 
we present in below. 

Because of 

1 g{t, 0)dx 



P[|x(t) - x(0)| > et] = J 



\x\>e 



vr (x - UQh{t, 0))2 + g^{t, 0) 



1 et + Uoh{t, 0) . et - Uoh{t, 0) 

= 1 — — arctan ; + arctan ; 

TT^ git,0) git,0) J 

we have 



,. „Jx t - X , _, 2 eX 

hm P\^—^ ^—^ > e = 1 arctan — 

i^oo ' t TT 0-2 



Therefore, ^ does not tend to zero even in probability. 

That means that sample paths of the process x(t) diverge to infinity faster 
than time t (up to dimensional constants). 



Remark 6: (i) We can generalize slightly the discussion and allow the 
parameter A to depend on time 

du{t) = -A(t)u(t)c/t + dB{t) 

where A(t) is a continuous function. Then, by integrating the above equation, 
we obtain 



u(t) = 6(t, s)u(s) + fb{t,T)dB{T) 



where 6(t, s) = exp[— / A(r)(ir]. By invoking our previous arguments it is 

s 

easy to find the probability distribution and transition function for the new 

process u(t) (and the new process of displacements). 

We can also consider an ra-dimensional situation, when B{t) is an R"- valued 
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Cauchy process, and the Langevin equation (1) is replaced by the following 
one 

du{t) = -Au{t)dt + dB{t) (32) 

where A denotes now n x n matrix with real coefficients. The existence of 
the solution for (32) in a general setting of ?i-valued processes, H being a 
real and separable Hilbert space, was established in |19 . 
(ii) Sample paths of the process x(t) are also cadlag. 



5.3 Interpreting w.{t) as a velocity variable for x(t): 
Limitations 

Finally, we shall discuss the relation between the process of velocities u(t) 
and the process of displacements x(t). In the Ornstein-Uhlenbeck case the 
process uouif) is continuous in the mean square, that is 

\im E[\uou{t + h) - uoc/(t)r] = 

h — ^0 

That follows from the continuity of its covariance function 
(ti, h) -^ E[uou{ti)uou{t2)] = aV^I*^-*2l 

Therefore, :s.ou{t) exists as a limit in mean square of the corresponding 
Riemann sums. Moreover, since sample paths of uou{t) are continuous, the 
integral exists also almost surely and they both coincide. Hence :x.ou(t) is 
not only mean square differentiable but is also differentiable in the sense of 
conventional mathematical analysis and its derivative is uou{t)- 

For the Cauchy process B{t) the situation is different. Because the mo- 
ments of x(t) do not exist so x(t) is not even continuous in mean square. 
Moreover, since its sample paths are not continuous, they have no deriva- 
tives either. 

However, x(t) has a velocity field u(t) in a probabilistic sense. 

Indeed, because u(t) is stochastically continuous so for /i — > there holds 

x(t + h) - x(t) 1 



h 

t 



— I u{T)dT — > u(t) 
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in probability, compare e.g. Eq. (17). 

That demonstrates that our naming (a priori) of u(t) the velocity ran- 
dom variable, can still be maintained to a certain extent, once we pass to 
the induced spatial variable and try to recover back a built-in information 
about the velocity process. This feature is slightly amusing and somewhat 
counterintuitive, since both the Cauchy noise-supported velocity process and 
the induced process of spatial displacements are discontinuos with probabil- 
ity 1. Anyway, the notion of velocity in the standard Ornstein-Uhlenbeck 
process has its own limitations as well: its nondifferentiability and thus the 
nonexistence of accelerations is not resolved but merely bypassed by invoking 
the white noise calculus. 

One should perhaps recall at this point that our theoretical framework 
reduces to a stochastic modelling of physical phenomena. That constitutes 
a fine-tuned approximation, in terms of stochastic processes, of a generically 
robust Reality. Surely, we have not attempted her genuine reproduction. 
The obtained description is much too detailed and thus necessarily involves 
a number of artefacts (the nonexistenece of accelerations in the standard 
Ornstein-Uhlenbeck process is one of them, an unbounded variation of the 
Wiener process sample paths just another). 

Acknowledgement: Both authors receive support from the KBN re- 
search grant No 2 P03B 086 16. 
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